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1. Introduction 

Colloidal suspensions of charged macroions — nanometers to micrometers in size and 
dispersed in a fluid by Brownian motion — are ubiquitous in nature and industry |lll2]- 
The remarkable thermal, optical, and dynamical properties of colloidal materials hinge 
on a delicate balance between competing interparticle interactions [3J. Common 
examples include aqueous paints, detergents, and clays: dispersions of latex particles, 
ionic surfactant micelles, and mineral platelets, respectively. Self-assembled crystals 
of synthetic, monodisperse, silica or polystyrene microspheres provide useful scaled-up 
models of atomic crystals and promise novel technologies, such as photonic band-gap 
materials [HE]. As predicted by the classic theory of Derjaguin, Landau, Verwey, 
and Overbeek (DLVO) [6l[7], repulsive electrostatic forces between charged colloids can 
stabilise a suspension against aggregation induced by van der Waals attractive forces [8]. 

Dispersed in a polar solvent, colloidal particles can acquire charge through 
dissociation of ionizable chemical groups at the surface. Electrostatic interactions are 
sensitive to the surface charges of the macroions and to the distribution of surrounding 
counterions. A macroion's bare (structural) charge depends on its surface chemistry 
(e.g., number and type of ionizable sites) and, in general, on the pH and salinity of the 
solution [9], Ho]. Dressed by an entourage of strongly attracted counterions, a highly 
charged macroion can act as though carrying a significantly reduced (renormalized) 
effective charge. 

The basic concepts of charge renormalization and effective charge were first 
introduced and widely explored some four decades ago in the context of polyelectrolyte 
solutions [m [12] . Similar ideas were subsequently applied to colloidal suspensions by 
Alexander et al. [13], who demonstrated that strong association of counterions can 
significantly renormalize spherical macroion charges. Numerous experimental studies of 
deionized aqueous suspensions of highly charged spherical latex particles [T^ fT5l fT6l fT7] . 
integral-equation calculations [18] , and simulation studies of the primitive model [HI [20l 
I2TI |22| [23] have since confirmed the effective charge as a physically important parameter 
in the one-component model of colloidal suspensions. 

The one-component model provides a practical approach to overcoming the severe 
challenges of extreme size and charge asymmetries in explicit molecular models of 
charge-stabilised colloidal suspensions, polyelectrolyte solutions, and many other soft 
materials [21]. The model is derived from the multi-component ion mixture by averaging 
over the degrees of freedom of the microions (counterions and salt ions). The surviving 
"pseudo-macroions" are governed by effective electrostatic interactions, screened by the 
implicitly modeled microions. 

This paper seeks to unite the concepts of dressed macroions and effective 
interactions in a coherent statistical mechanical framework that describes the association 
of counterions with macroions, the renormalization of the effective macroion charge, 
effective electrostatic interactions between macroions, and thermodynamic properties 
of deionized suspensions of highly charged colloids. Conceptually similar syntheses 
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have been proposed recently, based on Debye-Hiickel theory |2Sl ESI ESI ES] and 
on nonhnear Poisson-Boltzmann theory [30l [31] . The present theory is inspired by the 
elegant liquid-state approaches of Levin, Trizac, and coworkers [25l [261 EH [281 EH], but 
differs in several significant practical respects. 

The remainder of the paper is organized as follows. Sections [2] and [3] trace a 
path from the microscopic primitive model of charged colloids to an effective one- 
component model of dressed, charge-renormalized macroions. A simple criterion is 
adopted to differentiate between free and electrostatically bound counterions; physical 
approximations are developed for the free energies of the two counterion phases; and a 
variational method is prescribed for determining the renormalized effective charge and 
screening constant. Section [H demonstrates the practical implementation of the theory 
and compares predictions for the pressure of deionized suspensions with corresponding 
data from both simulations of the primitive model and experiment. Excellent agreement 
is obtained, over broad ranges of system parameters, with trivial computational effort. 
Finally, Sec. [5] closes with a summary and perspectives. 

2. Model 

Within the primitive model of charged colloids, the macroions are modeled as negatively 
charged hard spheres of monodisperse radius a and bare valence Zq (charge —Zqc), the 
microions as monovalent point charges, and the solvent as a dielectric continuum of 
uniform relative permittivity e. Polarization effects and image charges are ignored, 
assuming index- matching of macroions and solvent. The suspension may be either 
entirely confined to a closed volume at fixed salt concentration or in partial chemical 
(Donnan) equilibrium (e.g., via a semi-permeable membrane) with a microion reservoir, 
which fixes the microion chemical potentials. The reservoir is presumed to be a 1:1 
electrolyte solution with number density uq of monovalent salt ion pairs. 

In the presence of a sufficiently strong attractive potential, some fraction of 
counterions may remain closely associated with the macroions. By analogy with 
Oosawa's two-phase theory of polyelectrolyte solutions [12], a distinction then can 
be drawn between free and bound microion regions ("phases"). In contrast to 
rodlike polyelectrolytes, however, spherical colloidal macroions do not generate a 
Coulomb potential of sufficient range to overcome counterion entropy and condense the 
counterions. As a result, associated counterions remain only thermally (not physically 
or chemically) bound to the macroions. 

As Fig. [1] depicts, counterions localized within a spherical shell of thickness 6 (yet 
to be determined) are regarded as renormalizing the bare macroion valence. Colons are 
assumed to be completely expelled from the shell. The resulting "dressed" macroion is 
a composite object consisting of a bare macroion and its shell of bound counterions with 
an effective valence Z < Zq. Although the bare and effective valences are statistically 
fluctuating quantities, they are represented for present purposes by their average values. 
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Figure 1. Model of charged colloidal suspension: spherical macroions of radius a and 
point microions dispersed in a dielectric continuum. Strongly associated counterions in 
a spherical shell of thickness 6 renormalize the bare macroion valence Zo to an effective 
(lower) valence Z. 



The boundary between free and bound counterions is located at a distance from 
the macroion surface at which the electrostatic energy of counterion-macroion attraction 
is comparable to the average thermal energy per counterion. Denoting by 0(r) the 
electrostatic potential at distance r from a macroion centre, the association shell is 
defined via 



where C is an adjustable, dimensionless parameter, evidently of order unity. Counterions 
within the association shells {a < r < a + 6) are assumed to be trapped in the potential 
wells of the macroions, while more distant counterions have sufficient kinetic (thermal) 
energy to escape. This simple criterion for 6 justifies a Debye-Hiickel-like linear-screening 
approximation for the free counterions, which is exploited in the theory developed below. 

Previous studies have applied a thermal criterion similar to Eq. ([T]) to the 
electrostatic potential [321 [33l [31] or to the effective pair potential [35]. Alternative 
approaches to defining the association shell thickness are based on the structure of 
the counterions around a macroion or on the macroion configurations. Alexander et 
al. [13], for example, determined Z in a spherical cell model by matching the solutions 
of the nonlinear and linearized PB equations for the counterion density at the edge 
of the cell. The infiection point in the running effective charge of the macroion has 
been identified as another sensible boundary between free and bound counterions |36j . 
Yet another fruitful approach is to fit the effective one-component model (with a 
screened- Coulomb pair potential) to either the static structure factor measured in light- 
scattering experiments [Hj or the pair distribution function computed in simulations 
of the primitive model [T9l [201 [2ll [221 [23] . A useful comparison of various criteria for 
defining effective charges is provided in ref. |34j . 




(1) 
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3. Theory 

The theory proposed here for charge renormahzation and thermodynamics of colloidal 
suspensions requires modeling the electrostatic potential and the total free energy of 
the system. For this purpose, the most popular framework is the Poisson-Boltzmann 
(PB) theory [37], a mean-field approach that is especially well-suited to suspensions 
with monovalent microions, whose correlations usually can be justifiably neglected. 
Combining the exact Poisson equation for the potential with a Boltzmann approximation 
for the microion density profiles (as functions of position r), r;,-|-(r) = -n-o exp[=P'?/'(r)], the 
PB theory is based on the Poisson-Boltzmann equation 

= ft^osinhV^ , (2) 

where '0 = Peep is the reduced potential (vanishing in the reservoir), /? = l/^ksT) at 
temperature T, kq = y/STrXsno is the Debye screening constant, and Xb = (3e^/e is 
the Bjerrum length. Equation ([2]) must be solved together with appropriate boundary 
conditions: V'?/'|r=a = ZqXb/o? and Vip = either as r — oo — far from the 
macroions in a bulk suspension — or at r = i? in a symmetric cell of radius 
R. Neglecting macroion-macroion correlations, and all but asymptotically long-range 
microion-microion correlations, the corresponding Helmholtz free energy takes the form 

PF=J2 [ drn,(r)[ln(n,(r)A3) - 1] + -i- / dr , (3) 

where A is the microion thermal de Broglie wavelength and the two terms on the right 
side represent, respectively, the ideal-gas free energy due to microion entropy and the 
total electrostatic energy. 

At distances r for which !'?/'('") I ^ the right side of Eq. ([2]) may be approximated 
by an expansion about the reservoir potential [ip = 0) to linear order in ip. Anticipating 
applications to deionized suspensions of highly charged colloids, however, the microion 
densities are here expanded instead about the mean (Donnan) potential of the 
suspension ip [381 EH HO] : 

V^V = Ko[si^h'?/' -|- coshxp{ip — ^jj)] . (4) 

In a bulk suspension of macroions with bare valence Zq, the solution of Eq. (jlj), with 
boundary condition ilj'{r) — > as r ^ oo, yields the potential generated by a single bare 
macroion [IT] : 

ip{r) = -ZqXb — , r>a , (5) 

I + K,a r 

with the bare screening constant 



K = /toy cosh?/; = y47rA_B(n+ + n_) . (6) 

Here n± = no exp(=i='?/;) = N±/[V{1 — rj)] represent the mean number densities of 
microions in the free volume, i.e., the total volume V reduced by the fraction rj occupied 
by the macroion hard cores. Note that the screening constant k, depends implicitly on the 
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average density of macroions, since the global constraint of electroneutrality relates the 
numbers of macroions (A^m) and microions (A^±) in the suspension via Z^N^ = N^ — N_. 
Combining the linearized PB equation [Eq. (jl])] with a quadratic expansion of the ideal- 
gas free energy [Eq. ([3])] about the mean microion densities yields the corresponding 
linear-screening approximation for the one-body part of the free energy per macroion: 

,/^g..N«.AVll-f^-f (7) 

where Xi = Ni/Nm and the three terms on the right side account for, respectively, 
the microion entropy, macroion self-energy, and the Donnan potential energy of the 
microions. 

The PB theory proves to be formally equivalent to a class of effective-interaction 
theories that map the macroion-microion mixture onto a one-component model (OCM), 
by integrating over microion degrees of freedom in the partition function, and that 
neglect all but long-range microion correlations pll H2] . The effective- interaction 
approach has been variously formulated as density-functional [l3l HU HHl HHl W7\ , 
extended Debye-Hiickel [48], distribution-function [49l [50], and response [5ll [52l [53l 
[HU [55] theories — all fundamentally equivalent [21], aside from technical differences 
in the incorporation of excluded-volume effects. The OCM is governed by an effective 
Hamiltonian comprising a one-body volume energy and summations over pair and, in 
general, many-body effective interactions. 

Linearizing the PB equation about the mean potential [Eq. and the PB free 
energy about the mean microion densities, is completely equivalent in the OCM to 
linearizing the microion free energy about a reference system of neutral macroions 
embedded in an electroneutral microion plasma [12] and neglecting many-body effective 
interactions. Furthermore, the volume energy in the OCM turns out to be identical to 
the linearized PB free energy [Eq. (^^)]. A significant advantage of the OCM, however, 
is its natural incorporation of effective interactions between macroions. An additional 
contribution to the total free energy then comes from the effective (reduced) macroion- 
macroion pair potential [52], [53] 

(3v,s{r) = ZIXb — , r > 2a . (8) 

V 1 + Kay r 

Further progress requires uniting the OCM-based linear-screening theory of charged 
colloids with the charge renormalization model of Sec. [2J To this end, the total free 
energy is first separated, according to 

P -ffree ~l~ -^bound ~l~ Pfn i (9) 

into contributions from free and bound microions and from effective interactions between 
macroions, respectively. The linear-screening theory is then applied only to the free 
microions, whose free energy per macroion is approximated by [cf. Eq. ([7])] 

/3/free = E 5ain(n.A^) - 1] - f — - f , (10) 

~^ 2 1 + K[a + 0) 2 n+ + n_ 
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where x± = N±/Nm, N± are the numbers of free microions, n± = N±/[V{1 — fj)], and 
rj = ri{l + 6/a)^ is the effective volume fraction of the dressed macroions. Generahzing 
Eq. (I5l), the electrostatic potential around a dressed macroion of effective valence Z and 
effective radius a + 5 is given by 

„K(a+(5) 

= TTi(^ — 

with a renormalized screening constant 



/AnXBin+ + n-) . (12) 
The association shell thickness is now specified by combining Eqs. ([T]) and (fTTl) . yielding 

[l + k{a + 6)]{a + 6) ' ^ ' 

and solving self-consistently for 6 (given Z), noting that k depends implicitly on 6. 

The free energy of the bound counterions decomposes naturally into entropic and 
energetic contributions. The first contribution is the ideal-gas free energy of the bound 
counterions, given exactly by 

/3Fid = AttN^ dr r^nb{r) [in {n^{r)k') - l] , (14) 

where nfe(r) is the number density profile of bound counterions within the association 
shell and the integral covers the volume of the shell from inner radius a to outer radius 
a + 5. Although nft(r) could be obtained by solving the nonlinear PB equation [Eq. ([2])] 
(as in ref. [31]), the present study explores a simpler approximation, \\i[niy[r)h?) ~ 
In (nfcA'^), which yields 

/3Fid^Ar„(Zo-Z)[ln(n,A3)-l] , (15) 

where = {Zq — Z)/vs is the mean density of bound counterions in the association shell 
of volume Vs = (47r/3)[(a + 6)^ — a^]. The second contribution to the bound counterion 
free energy is the electrostatic energy F^i required to assemble the total charge of the 
dressed macroions — bare and bound charge — from infinity. An exact calculation 
would again require knowledge of the bound counterion density profile. Here we simply 
assume nf,{r) to be sharply peaked near r = a and take 

In the case of macroion charges below the renormalization threshold, Z = Zq and F^i is 
a trivial constant that is irrelevant for thermodynamics. At charges high enough that 
Z < Zq, however, Fgi becomes significant, since Z is state-dependent (as seen below). 
Combining Eqs. (fT5|) and (fT6|l . the bound-counterion free energy per macroion is here 
simply approximated by 



/3/bound ^ (Zo - Z) 



Inl^^A 
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For a given bare valence Zq, the effective valence Z is prescribed by minimizing 
with respect to Z the total microion free energy [sum of Eqs. ffTOl) and (fTTj) ] at fixed 
temperature and mean microion densities: 

y-Q^iffree + /bound) j ^ ^ ' ^^^^ 

The same variational prescription has been adopted by Levin et al. [25], [261 121! • It is 
easily shown that the minimization condition is equivalent to equating the chemical 
potentials of counterions in the free and bound phases, under the constraint that Z and 
5 are related by Eq. f|T3l) . The effective valence and corresponding shell thickness in 
turn determine the effective screening constant k via Eq. (fT2|) . 

Once the effective valence and screening constant are determined, the effective pair 
potential between dressed macroions follows as 

/ pka \ 2 -hr 

/35efr(r) = Z^Xb , r > 2{a + 5) , (19) 

y 1 + y r 

from which the macroion free energy Fm can be computed via liquid-state theory 
or computer simulation. (Note that the macroion radius is not renormalized in the 
prefactor of the effective pair potential in Eq. f|T9|) . since the association shells are 
penetrable.) For demonstration purposes, we implement a variational method [43l [55] 
based on first-order thermodynamic perturbation theory with a hard-sphere reference 
system [56j|. The macroion free energy per macroion is thus approximated as 

/m('^m,^±) = min(/Hs(n^,n±;c?) + 27m„ /" dr r'^g^s{r,nra]d)v^s{r,nrn,n±)\ , (20) 

where the effective hard-sphere diameter d is the variational parameter and /hs and gns 
are the excess free energy density and (radial) pair distribution function, respectively, 
of the HS fluid, computed here from the near-exact Carnahan- Starling and Verlet-Weis 
expressions [56]. Minimization of with respect to d generates a least upper bound to 
the free energy. It is important, in practice, to fix the renormalized system parameters 
[Z, 6, k) in this minimization and in all partial thermodynamic derivatives. 
The thermodynamic pressure finally can be calculated from 

P = ( ) " ^f^^^ + ^™ ' 

where / = F/Nm is the total Helmholtz free energy per macroion, Ns = is the 
number of salt ion pairs in the suspension, 

o ~ , ~ Z{n+-n^)k\B , . 

4[1 + K{a + 0)1'' 

is the (reduced) pressure generated by the free microions, and 

/?Pm = + ( ) (23) 

is the macroion pressure due to macroion entropy and effective pair interactions. Note 
that, since Z and S are implicitly held fixed in the partial derivatives, the bound 
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counterions make no contribution to the pressure. As an alternative to variational 
theory, computer simulation also can be used to determine the macroion pressure |57j . 

4. Results and Discussion 

To demonstrate its implementation, the charge renormalization theory is now applied 
to deionized suspensions of charged colloids and monovalent microions in an aqueous 
solvent at room temperature (A^ = 0.72 nm). As noted in Sec. [3l the theory involves 
a single free parameter, namely the dimensionless parameter C in Eq. ( |T8l) . which 
establishes the threshold for charge renormalization. To ensure that a counterion's 
average thermal energy does not exceed its binding potential, C must be 0{1). Lacking 
an independent, physical criterion, C must be regarded for the present as a fitting 
parameter. All results presented below were computed for C = 3, a value found to give 
satisfactory overall agreement with thermodynamic and structural data from primitive 
model simulations. In passing, we note that the thermal parameter C in the present 
theory is somewhat analogous to the adjustable cell radius parameter b in ref. [31], which 
combines PB cell and linear-screening theories. 

The key physical concepts of the charge renormalization theory are illustrated in 
Figs. [21111 For a sufficiently small bare valence, Eq. (fT3l) admits no real solution for 
nonzero thickness of the association shell. In this case, there are no bound counterions 
{5 = 0, Vs = 0) and the free energy is minimized hj Z = Zq (dashed line in Fig. [2|). 
At a threshold value of the bare valence, however, the shell emerges continuously and 
thickens rapidly with increasing Zq at fixed volume fraction and salt concentration (inset 
to Fig. [2|), while the free energy minimum shifts to Z < Zq (solid curve in Fig. [2|). The 
effective valence does not saturate, but continues to grow logarithmically with increasing 
Zq, in contrast to the behaviour expected and observed for polyelectrolytes [12] and to 
predictions for colloidal suspensions from PB cell-based theories [131 IMl EQ] and Debye- 
Hiickel-based theories |25l EHl EH EH [29] . 

As Figs. [3l and [H show, Z varies with volume fraction r] and salt concentration Cg of 
the suspension. The effective valence thus depends nontrivially on the thermodynamic 
state, exhibiting a pronounced minimum with respect to rj and increasing monotonically 
with Cs- Correspondingly, the association shell thins with increasing rj and Cg (Figs. [3] 
and [U insets), in such a manner, however, that the shells surrounding neighbouring 
macroions always remain separate and distinct. The absence of overlapping shells 
provides an internal consistency check on the theory. 

To test predictions of the theory for thermodynamic properties, the pressures 
of deionized suspensions, calculated from Eqs. ( fT2l) . ( fT3|) . and ( [2TI) -( !23l) . are directly 
compared with available data from simulations of the primitive model. Figures [5] and [6] 
show comparisons with the results of Linse [21] from extensive Monte Carlo simulations 
of salt-free suspensions with various bare valences and electrostatic coupling parameters 
r = A^/a. The unrenormalized linear-screening theory [52l [53] performs excellently 
for low-to-moderate couplings (F < 0.1779 for Zq = 40 in Fig. [5]), but breaks down at 
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Figure 2. Effective valence Z vs. bare valence Zq for a deionized suspension (cj ~ 0) 
of macroions of radius a = 50 nm and volume fraction rj = 0.1. Inset: counterion 
association shell emerges and thickens beyond threshold Zq. 
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Figure 3. Effective valence Z vs. volume fraction rj for a deionized suspension of 
macroions of radius a = 50 nm and bare valence Zq = 10"'. Inset: association shell 
(solid curve) thins with increasing 77, remaining thinner than nearest-neighbour surface 
separation in fee crystal (dashed curve). 



higher couphngs, characteristic of highly charged latex particles and ionic surfactant 
micelles. As illustrated in Figs. [5](b) and [6] (inset), charge renormalization becomes 
important for ZqT > 7, where the effective valence tends to be lower than the bare 
valence. The renormalized theory restores close agreement with simulation up to at least 
ZqT ~ 28 (r = 0.7115 in Fig. [5]). In practice, the excluded- volume correction to the 
microion densities in Eq. ( |T0|) . and the inclusion of the effective pair pressure — already 
important features of the unrenormalized theory [57] — are essential for consistent 
quantitative accuracy as the volume fraction becomes renormalized. Remarkably and 
intriguingly, the threshold for charge renormalization coincides with the onset of a 
spinodal phase instability, at low but non-zero salt concentrations, predicted by linear- 
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Figure 4. Effective valence Z vs. system salt concentration Cs for macroions of radius 
a = 50 nm. bare valence Zn — 10'*. and volume fraction n = 0.1. Inset: association 
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Figure 5. (a) Total reduced pressure Pp/ntot vs. macroion volume fraction rj, where 
ntot — {Zq + l)nm (total ion density), of salt-free suspensions with bare macroion 
valence Zq = 40 and electrostatic coupling constants (top to bottom) F = 0.0222, 
0.0445, 0.0889, 0.1779, 0.3558, 0.7115. Symbols: Monte Carlo simulations of the 
primitive model |21j (symbol sizes exceed error bars). Curves: variational theory with 
(solid) and without (dashed) charge renormalization. The double-ended arrow points 
to corresponding curves for T = 0.3558. The dashed curve for F = 0.7115 is ofF-scale, 
the pressure being negative, (b) Corresponding ratio of effective to bare macroion 
valence Z/Zq vs. 77 for F = 0.1779, 0.3558, 0.7115 (top to bottom). For F < 0.1, no 
renormalization is predicted {Z = Zq). 



screening theories [l3l H5l H8l [55] . This rather unusual prediction, however, has not 
yet been confirmed by primitive model simulations and the experimental situation is 
unresolved. 

Finally, the theory also can be tested against available experimental data. Figure 
[7] shows a comparison of predictions with osmotic pressure measurements of deionized, 
aqueous, charged colloidal crystals reported by Reus et al. [58] • Since bare (titratable) 
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Figure 6. Total reduced pressure (3p/ntot vs. electrostatic coupling constant F of 
salt-free suspensions with fixed volume fraction rj — 0.01 and bare macroion valence 
(top to bottom) Zo = 10,20,40,80. Open symbols: Monte Carlo simulations of the 
primitive model |21j (symbol sizes exceed error bars). Curves: charge- renormalized 
variational theory. Inset: Ratio of effective to bare macroion valence Z/Zq vs. F. 

charges are notoriously difficult to characterize in experiments, Zq is treated here as 
a fitting parameter, a value of Zq ~ 3000 giving a reasonable fit to the data. As 
seen in the inset to Fig. [Tj the effective charge in this case is substantially lower than 
the bare charge. It is important to emphasize, however, that only thermodynamic 
quantities have physical significance within the theory and that the theoretically defined 
variable Z does not necessarily correspond directly with any effective charge that may 
be determined experimentally, e.g., by light scattering, electrophoresis, or conductivity 
measurements. Moreover, direct comparisons between theory and experiment are 
subject to complication by charge regulation via chemical reactions at the macroion 
surface [9l [10], which may render even the bare charge dependent on thermodynamic 
state, e.g., pH and salinity. 
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Figure 7. Osmotic pressure (in Pa units) vs. volume fraction rj for a deionized 
suspension [cg — 0) of charged macroions of radius a — 51 nm and bare valence 
Zq =^ 3000. Curve: charge renormalization theory; symbols: experimental data [58j . 
Inset: effective valence Z vs. 77. 



5. Conclusions 

In summary, a new theory of charge renormahzation in charge-stabilised colloidal 
suspensions has been developed and implemented. The theory posits the existence of 
free and bound counterion phases and integrates a thermal criterion for distinguishing 
between the two phases with an effective-interaction theory based on a one-component 
model. Within the theory, bound counterions act to renormalize the effective valence 
of the dressed macroions, while free counterions screen the dressed macroions and 
make the dominant contribution to the pressure. A linear-screening approximation 
accurately describes monovalent free counterions, while the bound counterions are 
adequately described by a comparatively crude coarse-grained approximation for the 
bound-counterion density profile. 

Despite the conceptual and practical simplicity of the charge renormalization 
theory, predictions for the pressure closely agree with corresponding data from both 
primitive model simulations and an experiment, over ranges of macroion charges, volume 
fractions, and electrostatic coupling strengths, demonstrating the practical potential 
of the theory for modeling equilibrium thermodynamic properties. A preliminary 
simulation study [59] indicates that the theory also can accurately model structural 
properties, such as macroion-macroion pair distribution functions, although within a 
more limited range of electrostatic couplings. Future work will focus on refinements 
of the theory, further comparisons with experiment, and applications to the phase 
behaviour of deionized suspensions of highly charged macroions in bulk and in 
confinement [6OI EH, |62] . 
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